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Method for conjoint channel and direction of arrival estimation 

The present invention concerns in general terms a method of determining the 
spatial transfer function of a transmission channel by means of an array of antennae. 
More precisely, when the propagation of the signal between a transmitter and an array 
of reception of antennae takes place along a plurality of paths, the present invention 

5 makes it possible to estimate, for each path, the direction of arrival conjointly with the 
attenuation and phase rotation undergone by the signal along this path. The present 
invention applies more particularly to the estimation of the channel and the direction 
of arrival in mobile telephony. 

The propagation of a signal between a transmitter and receiver, in particular 

10 between a mobile telephone and a base station, can take place along multiple paths. 
Each path is characterised by a delay, a direction of arrival, an attenuation and a phase 
rotation. 

It is known that an array of antennae can be used to form a beam and/or cancel 
interference in one or more directions. The antenna processing consists of a weighting 

15 of the outputs of the different antennae but means of complex coefficients before 
summing, the coefficients being chosen so as to obtain the equivalent antenna diagram 
required. It is thus possible to form a beam in the direction of arrival of the useful 
signal whilst placing zeros in the directions of arrival of the interfering signals. The 
majority of the beam formation techniques however require prior knowledge of the 

20 direction of arrival of the signal. Beam formation has been applied to mobile 
telephony, notably to direct a reception beam from a base station to a mobile station 
(uplink). The base station is then equipped with an adaptive antenna (referred to as an 
"intelligent antenna") capable of pointing in the direction of a propagation path 
issuing from a mobile terminal. 

25 In addition to the direction of arrival of the useful signal, it is often necessary to 

know the characteristics of the transmission channel. Thus, when a mobile telephony 
system is functioning according to a code division multiple access (CDMA) mode and 
is exploiting diversity in reception by virtue of a rake receiver, the complex coefficient 
used in the different branches of the receiver are obtained from the attenuation and 

30 phase rotation undergone along the different paths. 

It is known that an adaptive antenna can be used, exploiting the diversity of 
paths. For this purpose, the systems of the state of the art use one rake receiver per 
antenna and effect the beam formation from signals issuing from the different 
receivers. It can also be envisaged pointing several beams in the directions of arrival 
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of the different paths and then combining the different paths in a rake receiver. In 
both cases the direction of arrival of the useful signal and the characteristics of the 
transmission channel are estimated independently. This approach has the double 
drawback of requiring two estimation steps and not using all the information available 
5 at the second estimation step. 

The aim of the present invention is to remedy the aforementioned drawbacks 
and to propose a conjoint method of estimating the direction of arrival and the 
characteristics of the transmission channel. More particularly, when the transmission 
channel has several propagation paths, the invention makes it possible to estimate the 

10 direction of arrival of each path and the complex coefficient (attenuation and phase 
rotation) characterising each path. 

The invention is defined by a method of estimating the channel and the direction 
of arrival of a signal transmitted by a transmitter and received by an array of antennae 
after being propagated along at least one path, the method comprising, for each path, a 

15 first step of estimating phase differences (£ f ) in the signals received by the different 
antennae in the array, a second step of estimating the angle of arrival (6) of the signal 
as well as the phase rotation (v) undergone by the signal along the said path from the 
said phase differences and a third step of estimating the attenuation (a) undergone by 
the signal along the said path from the estimated values (v , 0) of the phase rotation 

20 and the angle of arrival. 

Advantageously, the first step comprises, for each signal received (xt), the 
minimisation of a first distance between a plurality of values of the said signal ( xe(t) ) 
taken at a plurality of times (r) and the values taken at the same times of a pilot signal 
(b(t)) which underwent the phase shift and the attenuation to be estimated. 

25 According to one characteristic of the invention, the second step comprises a 

removal of ambiguity of the phase differences, the removal of ambiguity operating 
step by step from one antenna to the next starting from a reference antenna for the 
array. 

Advantageously, the removal of ambiguity uses an affme relationship between 
30 the phase differences and the rank of the antenna in the array. 

According to another characteristic of the invention, the phase rotation (v) and 
the angle of arrival (0) are estimated from a linear regression on the estimated values 
of the phase differences. 
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Advantageously, the phase rotation (v) and the angle of arrival (0) are estimated 

by minimising a second distance J(v,cp) = J](v + (^-l).(p-t, J where |r is the 

r=i 

estimated value of the phase shift of the signal received by the antenna of rank I , 
<p =2n.cos(6).dlX where d is the pitch of the array, X the wavelength of the signal and 
5 L the number of antennae in the array. 

According to another characteristic of the invention, a new estimation of 
the phase differences is carried out using the estimated value of the phase rotation (v) 
and that of the angle of arrival (#). 

Advantageously, the attenuation (a) is estimated by minimising a third distance 
10 between a plurality of values of the said signal (xt(t) ) taken at a plurality of times (r) 
and the values taken at the same times of a shifted pilot signal (b(t)) of the phase shift 
values (|r) newly estimated and which have undergone the attenuation (a) to be 
estimated, the distance being calculated on all the antennae in the array. 

The invention is also defined by a signal reception device comprising an array of 
15 antennae and estimation means adapted to implement the steps of the estimation 
method disclosed above. 

According to one characteristic of the invention, the device comprises, at the 
output of each antenna, a plurality of filters adapted to the different signal propagation 
paths and a plurality of estimation means, each estimation means being associated 
20 with a path (i) and receiving the outputs of the filters adapted to the corresponding 
path. 

According to another characteristic of the invention, the device comprises 
channel formation means, each channel formation means being associated with a path 
(i) and receiving, from the estimation means associated with the said path, the 

25 estimation 6 \ of the angle of arrival of the said path in order to form a channel in the 
corresponding direction. 

In addition, each channel formation means can also receive estimation means 
other than the associated estimation means, the estimations (6 V , iVi) of the angles of 
arrival of the other paths in order to place zeros in the corresponding directions. 

30 Advantageously, the device comprises a plurality of complex multiplication 

means, each complex multiplication means being associated with a path (i) and 
multiplying the output of a channel formation means by a complex coefficient 
di.Q~ J( ' where v \ and a \ are the estimated values of the phase rotation and of the 
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coefficient of attenuation supplied by the estimation means associated with the said 
path. 

The characteristics of the invention mentioned above, as well as others, will 
emerge more clearly from a reading of the description given in relation to the 
5 accompanying figures, amongst which: 

Fig. 1 depicts schematically a multipath propagation between a source and an 
array of antennae; 

Fig. 2 depicts a diagram illustrating an ambiguity removal procedure useful to 
1 0 one embodiment of the invention; 

Fig. 3 depicts a reception device according to one embodiment of the invention. 

It is assumed, as illustrated in Fig. 1, that a signal is transmitted by a source and 
received by an array of antennae. The array illustrated is uniform and linear but, 

15 naturally, other embodiments can be envisaged, notably sectorial arrays having 
elementary antennae arranged in an arc of a circle. The signal propagates according to 
a plurality p of paths, each path i=l..p arriving at the array at an angle of incidence 9j. 
The set of paths constitutes the transmission channel between the transmitting source 
and the receiving array. If several sources j=l..J each transmit a signal Sj, there will 

20 be J distinct transmission channels, each channel consisting of a plurality of paths. 

In the context of a mobile telecommunication system, the transmission channels 
in question will correspond to uplinks between mobile terminals, hereinafter referred 
to as users, and the base station. When several uplinks are simultaneously established 
with the same base station, a given transmission channel may be disturbed by the 

25 signals from the other users (multi-user interference). In addition, the signal of a 
given user being propagated on different paths arrives at the array after different 
propagation times. The transmission channels generally being dispersive in time 
because of their frequency selectivity, their pulse responses have a not insignificant 
duration, which creates interference between the signals which have come from the 

30 same user but have followed distinct paths (multi-path interference). In certain 
systems, notably those using direct sequence CDMA (DS-CDMA) access, each path 
of each user comes to be isolated to a certain extent. To do this, the good properties of 
autocorrelation and intercorrelation of the spreading sequences (spreading codes) are 
exploited by filtering the signal received by a bank of suitable filters. The multi-user 
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and multi-path interferences are however not completely eliminated and residual 
interference remains. If the users are sufficiently numerous, the limit central theorem 
can be applied and the residual interference seen an additive white Gaussian noise. 

The invention presupposes that the receiver can, at least roughly, separate, at 
5 each antenna, the different paths of a signal issuing from a user. This will be the case 
for a DS-CDMA system, if a battery of filters adapted to these different paths is used 
at the output of each antenna. This will still be the case if the transmission channel 
between the user and receiver consists essentially of a single path. 

Each path i, i=l..p is in reality the combination of a large number of elementary 
10 paths which have undergone reflections and diffusions produced by the environment 
of the receiver and transmitter. With each of these elementary paths there are 
associated a phase rotation and an attenuation which follow uniform laws. According 
to the limit central theorem, the signal s,(t) coming from the path i can then be 
written: 

15 

5,(0 = a,{t).e JV ' t0 .6,(0 + «,(0 (1) 

where a,{t) is the attenuation associated with the path i obeying a Rayleigh law, 
v.(0 is the phase rotation, uniformly distributed on ]->r;?r], 6, (0 is the information 

20 sent and n,(0 is an additive white Gaussian noise including the residual interference. 
The channel is variant over time. One of the parameters making it possible to quantify 
this variation is the coherence time (A0 C , which guarantees that two responses of the 
channel at two times separated by more than (At) c are non-correlated. If the antennae 
employed in the array are not omni-directional, the coefficient «,.(0 also includes the 

25 antenna gain. 

The transmission channel between a user and the receiver will now be 
considered. It will be assumed that the receiver can isolate a path i of the channel, 
either because the channel includes only one path, or because a filtering adapted to the 
path i is used. A single path being considered, it is possible to omit the index i in 
30 order to simplify the notations. It will also be assumed that the signal arrives at a 
uniform linear array of L antennae at an angle of incidence 6(0 . 

Denoting as x e the signal received at the antenna t, i = (1, . . . , L) , it is possible 

to write: 



x e (t) = a(t).e J4 ' 0) .b(t) + n ( (t) 



(2) 



with 

^(0=K0+(^-i)-?<0=K0+2^-i)cos(6(0) . (3) 

where d is the pitch of the network, X is the wavelength, b(t) is the symbol sent 
by the user, n t (t) is the additive white Gaussian noise at the antenna t including the 
thermal noise, the multi-path interference and the multi-user interference. 

The response of the channel is thus modelled by an attenuation coefficient a(t) 
and a phase rotation v(t) . The total phase shift £,(0 results on the one hand from the 
phase rotation v(/) caused by the channel and on the other hand from the difference 
in operation between the antenna I and the reference antenna 1 at the angle of arrival 

The invention proposes a conjoint estimation of v(t) , 6(f), a(t) from the 
signals x t (t) and a priori information on the transmitted symbol b(t). The symbol b(t) 
can be a pilot symbol perfectly known to the receiver or a symbol issuing from a 
decoding operation prior to soft or hard decisions and having a sufficient degree of 
reliability. 

v(0 , 6(0 , ct(t) will be assumed constant over an observation window of length 
T chosen so as to be small compared with the coherence time (At) c . The length T 
will be chosen so that, the higher the speed of the mobile, the shorter the length, and 
chosen so that, the lower the signal to noise ratio, the longer the length. It can change 
over time according to the transmission conditions. 

0 , a and v are sought such that the quadratic distance D is at a minimum. 

The resolution of such a problem requires the resolution of a complex non-linear 
system. The problem is resolved by dividing it into three steps. 

Step 1: 4 t is estimated by minimising the quadratic distance D t at each 
antenna: 
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D t =^\x t {t)-a.e^.b{t^ (5) 
Step 2: 6 and v are estimated from a linear regression on the estimated values 

i t - 

5 

Step 3: from the estimated values v and 6, a is estimated finally by 
minimising the quadratic distance 

In a first step it is sought to estimate the phase shift of the signal at each 
10 antenna, that is to say the value £ t of & which minimises D t : 

= £ ( Xl (o - a** . w)- (*; (o - «■«-*« •** (0) 

= Z( M')f + « 2 -|W| 2 - 2.Rc(r £ (/)^-* i>*(0) ) 

(6) 

If the real part of equation (6) is taken, there is obtained: 
x ( {i).a.e*> .b" (0 = a.(xf (/) + ./.^ (0) V (0 " (0)(cos(£) - y.sin(^ )) 
15 =«.([xf (0-^(0 + ^(0-^ «]+ Ax' t W*(t)-xH0.b I (t)] ) 

.(cos(£)-y.sin(£)) 

(7) 

and therefore 

Mx t {t).a.e j s> b * (0 ) = a ( [ x f (0>ft * (,) + , 4 ' (0j.cos(^) 

+ [xUt).b\t)-x?(t).b'(t)].sm&)) 

(8) 

20 By using equations (6) and (8) it is possible to calculate the derivative of the 

distance D ( according to . There is obtained: 
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dD t = 2 ^ d\Re{x ( (t).a.e- J? '.b\t))] 

= - 2 «-S ( - (0 + A it).b ] (0]-sin(^) + (0 V (/) - Jcf (/).*' (0}cos(^ ) ) 

(9) 

The value £, t of is sought which cancels out the derivative of equation (9). 
4 t must therefore satisfy the following equation: 

/=! /=1 

(10) 

In order to simplify the writings, let the following sums be defined: 
(=1 

By replacing the above sums in equation (10), the estimated value of the phase 
10 shift £ e is obtained: 

01) 

In the second step, 9 and v are estimated from the estimated values % t of the 
phase shifts % t . The calculation of the estimated values 6 and v will be made in the 
15 simple case of a uniform linear array but can extend to any array. 

On the first antenna ( I — 1 ) taken as the reference antenna, the phase rotation 
does not depend on the direction of arrival: = v . It is therefore possible to 
provisionally choose v equal to . Equation (10) nevertheless has two distinct 
solutions on ~\-7t\7t\, which correspond to a maximum and minimum of the function 
20 D t {£ ( ) and different by an angle n . These two estimated values are: 

| 1 0 e^-|;|jand| 1 1 =| 1 °+^ (12) 
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The minimum corresponds to the value for which the second derivative 5 D ' with 
£ = 1 is positive. 

^ = 2als] + S< ).cos(<f f ) + (s\ - S] ).sin(£,)] (13) 

|, is that of the two values supplied by equation (10) which corresponds to a. positive 
value of formula (13) for I = 1 . Consider now the other antennae {i > 1 ). The value 
of £ f is related to 0 and v by the following equalities: 

4 t =v + p t = v + (£-l).<p = v + 2xj-(£-l)cos0 (14) 
Equation (10) has several solutions. Let £° be the unique solution such that 
<f° e ^ ; yj . The other solutions are = |° + k.n with & a relative integer. 
According to equation (14) there must always be 

\\ e y2n±(£-l) + v;2n*(t-l) + v^ (15) 

It must therefore always be necessary to seek all the possible solutions which satisfy 
conditions (14) and (15) replacing v by its estimated value already available If 
d I X = 1 12 , the length of the segment of equation (15) is 2(i - \)ti . There therefore 
exists at most 2{l — l) possible values of for d 1 1 < 1 / 2 . It is then necessary to 
verify, for each of these solutions, if it is indeed a minimum and therefore if the 
second derivative supplied by formula (13) is indeed positive for the value % t = • 
By virtue of the properties of periodicity of the cosine and sine functions, it can be 
shown that, if |° is a minimum, all the |* values with k even are minima and all the 
|* values with k odd are maxima. Likewise, if £° is a maximum, all the values 
with k even are maxima and all the values with k odd are minima. It is therefore 
deduced from this that, amongst the l{l - 1) possible values of |* , only (t - 1) values 
are minima. 

It therefore remains to choose one minimum amongst the different minima obtained . 
The ambiguity removal procedure is illustrated in Fig. 2. For £ = 2 and d I X < 1/2 , 
which is always the case in practice, |* € \- it + v ; n + v] and there therefore exists 
only one minimum | 2 . This minimum makes it possible to make a first estimation 
(p m of <p: <p {2) =£ 2 -^- The ambiguity on the estimated values |* increases 
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linearly with the rank £ of the antenna. However, the affine relationship between 
and q> and a condition of continuity of (p from one antenna to the following make it 
possible to remove this ambiguity. More precisely, for £ > 2 , amongst all the possible 
minima the one is chosen which leads to the smallest variation in the estimation of <p 
from the antenna £ - 1 to the antenna £ , that is to say: 

| ( =|* with Jt= argmin ( |£(*.o -£«-d| ) : (16) 

k such that |? is a minimum 

where ^""fff 

where the choice of |, expressed by equation (16) represents the continuity of <p 
from one antenna to the next. The ambiguity is thus removed from step to step, from 
one antenna to the following one, starting from the reference antenna ^=1 . 

There is thus available a set of L estimated values £ t for £ = {\,...,L\ from 
which the estimated values v and y are then obtained by linear regression. More 
precisely, the values v and q> are sought which minimise the standard deviation J : 

J(v,<p) = f,{v + (£-l).<p-i} ( 17 > 
These two values are obtained by cancelling out the partial derivatives of the function 
J: 

— (v,#) = 2.y(v + -&)=<> <=* v = % -$1^ O 8 ) 

dv 7^ 



with f = and *-l = ■}-;£</- 1) 



20 ^(v,<p) = 2.Y(£-l){<> + (t-\).<p-i)=0 o {l-\f.$ = {£-\\?;-vt-\ 

(19) 

with ( 



with (T^iK = t2> ~ ^ and (£ " 1)2 = 7 2 ( ^ " 1)2 

By replacing equation (18) in equation (19), there is obtained : 



(£ — \\ P - £ £-\ 

-l) 2 .^ = (F^-(£ -<p.£-l)l-\ « <P = -— - I __ 2 - (20) 



25 
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By means of equation (18), it is then possible to calculate v . From the value of <p , it 

is also possible to deduce 6 : 

0 = arccosf— f—W<) (21) 
^ 7.7m I a) 

At the end of the second step estimated values v and 6 are available. In the 
third step it is sought to evaluate the attenuation coefficient a . To do this, first of all 
a new estimation % t is carried out of the phase differences at each antenna t using v 
and (f> , that is to say: 

f f= v + ^-l).^ ( 22 ) 
Finally, the estimated value a of the attenuation a of the channel is that which 
minimises D, knowing that & =4„V£. According to equation (4), the value a, 
dD 

which cancels out the derivative — — , satisfies: 
da 

v, ±±Re{x e «).e-*.b\t)) 

^(«) = WW!' -Re(*,(0.e"*' .6*(0)] = 0 ~ « = M t= ' , 

da tfM Xj;j6(/)| 2 

(23) 

The estimation of 9, v and a can be carried out dynamically using a sliding analysis 
window. Pursuing the angle of arrival then enables the beam former to continuously 
follow the mobile. Likewise, pursuit of the parameters v and a makes it possible to 
ensure an adaptive combination of the different branches in a rake receiver. 

Up to now a single path of the transmission channel has been considered. If 
the channel has several paths i=l..p, the conjoint estimation method must be applied to 
each path i. The channel can then be characterised by a set of triplets of estimated 
values (6\,v\,& s), i=l..p, each triplet itself characterising a path. 

Fig. 3 illustrates a reception device according to one embodiment of the 
invention. The device comprises an array of antennae (300,) ...(300 L ). Each antenna 
(300^) is connected to a battery of suitable filters (310s), i=l..p, each filter (3100 
being adapted to a path i. The outputs of the different filters (310s) relating to the 
same path i are connected to a conjoint estimation module (3300 implementing the 
estimation method described above and supplying a triplet of estimated values {9 h vx, 
a 0. The estimation is effected in parallel for all the paths. The output signals 
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relating to the same path i=l..P are also directed to channel formation means 
(320i,..,320 P ). The estimated value Q \ supplied by the conjoint estimation module to 
the channel former (320,) enables the latter to point its reception beam in the direction 
of arrival of the path i. In a variant embodiment which is not shown, each channel 
5 former (320j) also receives the estimated values 6\- supplied by the other estimators 
(330iO, iVi, of the paths i' of the same user. This enables the channel former (320j) to 
place zeros in the reception diagram for the directions concerned (0v, iVi) so as best 
to separate the signals received from different paths. The estimated values v ■, and d j 
are transmitted to a complex multiplier (3400, which multiplies the output signal of 
10 the channel former (3200 by the complex coefficient ai^' jyi , the conjugate of the 
complex coefficient of the channel. The outputs of the P multipliers are then summed 
in a summer (350) in order to supply an output of the MRC (Maximum Ratio 
Combining) type. 

This reception device makes it possible to effect both a formation of channels in 
15 the directions of the different paths and a combination of the MRC type as in a rake, 
filter. It is clear to a person skilled in the art that the operations represented by the 
different modules in Fig. 3 can be effected by a processor programmed for this 
purpose or by a plurality of dedicated processors. 

20 Although the invention has been essentially described in the context of the 

reception of a signal transmitted by a mobile terminal and received by the adaptive 
antenna of a base station, it is clear to a person skilled in the art that the invention is 
not limited to this example application but extends to any type of application requiring 
the estimation of the direction of arrival of a useful signal and the characteristics of a 

25 transmission channel. 



